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Late time dynamics of scalar perturbations outside black holes. 
II. Schwarzschild geometry 
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We apply a new analytic scheme, developed in a preceding paper, in order to calculate the late time 
behavior of scalar test fields evolving outside a Schwarzschild black hole. The pattern of the late 
time decay at future null infinity is found to be the same as in the shell toy-model studied in the 
preceding paper. A simple late time expansion of the scalar field is then used, relying on the results 
at null infinity, to construct a complete picture of the late time wave behavior anywhere outside the 
black hole. This reproduces the well known power-law tails at time-like infinity and along the event 
horizon. The main motivation for the introduction of the new approach arises from its applicability 
to rotating black holes, as shall be discussed in a forthcoming paper. 
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I. INTRODUCTION 

It is well established that the gravitational field of a 
generically forming black holes relaxes at late time to 
a “no-hair” stationary Kerr-Newman geometry. It was 
first demonstrated by Price j|], regarding gravitational 
and electro-magnetic perturbations of the Schwarzschild 
black hole (SBH) exterior, that the fields die off at late 
time with an inverse power-law tail. For a spherical- 
harmonic wave mode of multipole number I, it was shown 
that a i-( 2 U 2 -i-p) (tecay tail {t being the Schwarzschild 
time coordinate) will be detected by a static observer 
outside the black hole, with p = 1 if an initially compact 
perturbation is considered, or p = 0 in case that a static 
field existed outside the central object before the onset 
of collapse. 

These results were later confirmed using several differ¬ 
ent techniques, both analytic and numerical M , and 
were generalized to other spherically-symmetric space- 
times |H,||,0. The application of perturbative (linear) 
approaches is encouraged by numerical analysis of the 
fully non-linear dynamics of the fields |Q,H, which indi¬ 
cates virtually the same late time pattern of decay, as for 
the minimally-coupled (linear) fields. 

Power-law decay tails are exhibited by fields at late 
time, because in a curved spacetime waves do not prop¬ 
agate merely along light cones, even when the fields are 
massless. Rather, the waves spread inside the light cones 
due to scattering off spacetime curvature. As suggested 
by previous studies, the late time behavior of these waves 
is characteristic of merely the large distance structure of 
spacetime. This implies that the phenomenon of late 
time tails may not necessarily be restricted to the ex¬ 
teriors of black holes. For example, late time tails are 
found to form during the purely spherical collapse of a 
self-gravitating minimally coupled scalar field, even when 
the collapse fails to create a black hole § . Conversely, no 
power law tails are detected in the non-asymptotically- 


flat geometries of Schwarzschild-de Sitter and Reissner- 
Nordsrom-de Sitter black holes ^ (instead, the field is 
found to die off exponentially at late time in theses cases). 

All previously mentioned studies were benefited from 
the simplicity of spherical symmetry. Yet, an astrophysi- 
cally realistic model should clearly employ a rotating cen¬ 
tral object. Thus, apparently the most tempting gener¬ 
alization of the analysis involves the inclusion of angular 
momentum in the background geometry. A first progress 
in this direction has been achieved recently with the in¬ 
troduction of a full (1-1-2 dimensions) numerical analysis 
of wave dynamics in Kerr spacetime, by Krivan et al. 
| p^ , pl| . So far, however, no analytic scheme has been 
proposed for the investigation of wave dynamics in Kerr. 

In a preceding paper (to be referred to as paper I) 
we introduced an analytic technique for the study of late 
time behavior of fields in asymptotically-flat spacetimes. 
The prime motivation for the introduction of the new 
scheme was its applicability to rotating black holes. To 
examine the essential features of the proposed calcula¬ 
tion scheme, we applied it in paper I to study the simple 
toy-model of a scalar field evolving outside a spherically- 
symmetric thin shell of matter. In that case, the new 
technique, based on what we called “t/ie iterative expan¬ 
sion'' , allowed a simple and rigorous derivation of the late 
time waves-form at null infinity. In the present paper we 
apply a variant of the iterative scheme in order to analyze 
the evolution of scalar waves on the background of the 
complete Schwarzschild geometry. Again, this method 
will enable the analytic calculation of the late time be¬ 
havior at null infinity. We shall than show how, relying 
on the results at null infinity, it becomes rather simple 
to construct a complete picture of the late time decay 
anywhere outside the black hole, in particular along the 
event horizon. 

There are several reasons why we think it is worthwhile 
to first analyze the already well-studied case of a SBH, 
rather than directly focus on the more interesting case 
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of the Kerr black hole. First, this will enable us to test 
our scheme against the well-established results available 
in the Schwarzschild case. Secondly, many parts of the 
analysis in Schwarzschild shall later be directly employed 
when analyzing scalar waves in Kerr |^,|^. Finally, the 
analysis in Schwarzschild will appear to be valuable on 
its own right, providing, in some respects, a more com¬ 
plete picture of the late-time wave behavior than already 
available. 

In the shell model, spacetime is flat at small distances 
(inside the shell). For that reason, the complete internal 
geometry could be exactly accounted for by merely the 
“Minkowski-like” first component of the iterative expan¬ 
sion (denoted in paper I by dto). (we remind that the 
terms 'kAr>i of the iterative expansion were describing 
deviations from flat geometry, namely curvature effects 
outside the shell.) The complete Schwarzschild manifold, 
however, does not share this convenient property, as in 
this case spacetime is highly curved at small distances. 
This will enforce us to choose for another “basis” poten¬ 
tial for the iterative scheme at small r (other than the 
purely-centrifugal potential Vq chosen in the framework 
of the shell model), and will thus somewhat complicate 
the technical details of the analysis. Nevertheless, the 
basic calculation scheme, as well as the results at null 
infinity, shall remain essentially the same as in the shell 
model. 


arrangement of this paper 


This paper is arranged as follows. In sec. || we give 
a mathematical formulation of the wave evolution prob¬ 
lem in Schwarzschild as a characteristic two-dimensional 


initial-value problem. In sec. Ill we introduce the itera¬ 


tive scheme to be used to allow an analytic treatment of 
the mathematical problem. We apply the iterative cal¬ 
culation scheme is sections 0 through |VII| , obtaining an 
expression for the waves-form at late time at null infin¬ 
ity. Then, in sec. EM a simple technique is applied to 
obtain the late time behavior of the scalar field at any 
constant r (including along the event horizon). Sec. IX 
summarizes the results and discusses possible extensions 
of the analysis. 


II. THE INITIAL VALUE PROBLEM 

We consider the evolution of initial data, representing 
a generic pulse of massless scalar radiation, on a fixed 
SBH background. The scalar field is assumed to satisfy 
the (minimally-coupled) Klein-Gordon equation 

= 0 , ( 1 ) 

where $ represents the scalar wave. The structure of 
spacetime affects the evolution of the scalar field through 


the covariant derivatives, denoted in Eq. (|^) by semi¬ 
colons. 

Decomposing the field into spherical harmonics, 

OO I 

EE (2) 

m— — l 


we obtain an independent equation for each of the com¬ 
ponents r), 




, 2(r-M) . , ^(^ + 1)^, 


= 0 . 

( 3 ) 


Here t, r, 6 and ip are the standard Schwarzschild coordi¬ 
nates, M is the mass of the black hole, /(r) = (1—2M/r), 
and I is the multipole number of the mode under consid¬ 
eration. 

A more convenient form for the wave equation 
may be obtained in terms of a new wave function 
r) = r). To that end we introduce the double¬ 

null (Eddington-Finkelstein) coordinates v = t + and 
u = t — r*, where 


r.=r + 2Mln(^l^). (4) 

The ’’tortoise” coordinate r* varies monotonically from 
—OO (the event horizon) to -l-oo (space-like infinity). 

The wave equation now reads 

¥^, + V^{r)¥ =0, (5) 


in which 

V\r) 


= - I- 


2M 


l{l + 1) 2M 


( 6 ) 


is an effective potential, accounting for both centrifugal 
and curvature effects. This effective potential is sketched 
in figure || as a function of r* for the sample values I = 
0 , 1 , 2 . 



r*/2M 


FIG. 1. Effective potential for scalar waves in Schwarzschild 
spacetime. With r+ defined as in Eq. (Q), we have r* = 0 corre¬ 
sponding to r ~ 2.56M. 
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Note the following features of the effective potential 
(valid for all values of /), which play an important rule 
in our analysis. V(r) is localized (in a sense apparent in 
figure forming an effective potential barrier for the 
waves. At large distance, V(r) is dominated by the cen¬ 
trifugal potential (reflecting asymptotic flatness), with 
curvature-induced deviations which die off as ^ r~^. At 
small r* values, V(r) dies off exponentially in r*/M to¬ 
wards the event horizon, making the potential effectively 
zero inside the potential barrier. Evidently, the late time 
behavior at null infinity is affected mostly by the shape of 
the potential at large distance. Conversely, the early evo¬ 
lution (e. g. the quasinormal ringing stage), is strongly 
related to the fine details of the potential shape at small 
r values. 

Since each of the spherical harmonics modes evolves 
separately, we henceforth discuss the evolution of a single 
mode of arbitrary multipole number 1. The superscript 
‘Z’ (denoting Z-dependence) will usually be suppressed for 
brevity. 

The initial data for the evolution problem shall be 
specified on two characteristic (null) surfaces outside the 
event horizon, as sketched in the Penrose diagram of fig¬ 
ure ||. We will first consider initial data in the form of 
some compact outgoing pulse, specified on the ingoing 
null surface v = 00, 

/ '^{u = uo) = 0 , , 

I 4-(^; = 0 ) = r(w) ’ 

where r(M) is some function of a compact support (“the 
pulse”) between retarded times u = uq and u = ui. As 
demonstrated in paper I for the shell model, the case of 
static initial field can be later inferred in a simple way 
from the result regarding a compact pulse. 


i+ 



*As long as uo remains a free parameter in the analysis, 
the specific choice v = 0 causes no loss of generality, since 
spacetime is time-translation invariant. 


FIG. 2. The set-up of initial data. Shown is the Penrose dia¬ 
gram representing the external Schwarzschild geometry. The dark 
feature (artificially shown as if extended to u > 0 values) represents 
the amplitude of some compact support initial function r(tt) on the 
ray v = 0. The initial problem for the scalar field is well posed in 
region S (the shadowed area). 

The evolution equation (H), supplemented by the ini¬ 
tial conditions (|^, establishes a well posed characteris¬ 
tic initial value problem for the scalar field anywhere in 
the domain S outside the event horizon (see figure 2). 
Since, manifestly, this problem poses no mathematical 
irregularities, existence and uniqueness of a solution are 
guaranteed by fundamental mathematical theory (see, for 
example, [p^). 


III. THE ITERATIVE EXPANSION 

To define the iterative expansion to be applied in the 
complete SBH model, we first introduce a new parameter. 
To > 0, its value chosen so that r(r* = rg) is of order 
> 2M [say, r(r^, = rg) = 3M]. We then define 

{ 0 , r* < rg 

K'+i) „ > „ (8) 

and 


SV(r) = V-Vo, (9) 

in which V(r) is the Schwarzschild effective potential 
given by Eq. The potential Vg(r,) is so defined to 
account for the fact that the actual effective potential is 
exponentially small at small r* values. With this defi¬ 
nition, the function Vg(r*) approximates the form of the 
actual effective potential V (r*) at both the very large and 
the very small values of r* (The deviations, described by 
SV, become significant only at intermediate distances, 
see figure |l|) . 

Following the same procedure as in studying the shell 
model, we define the iterative expansion by decomposing 
the scalar wave dl into an infinite sum, 

CXD 

^ ( 10 ) 

N=0 

in which the components Ikare defined in a recursive 
way by the hierarchy of equations 

T , T f 0 , A = 0 

+ “ I _(5V)«'Ar_l ,A > 0 , 

supplemented by the initial data 

«'Ar(M = ug)=0 (VA > 0), (12a) 


'^n{v 



, A = 0 
,A > 0 ■ 


(12b) 
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Formal summation over N recovers the “complete” initial 
value problem for the scalar wave 

It was indicated in paper I that in the analogous shell 
model the iterative sum seems to converge rather effi¬ 
ciently at late time to the actual field at null infinity, 
provided that the initial pulse is specified at large dis¬ 
tance. In that case, it was suggested both numerically 
and analytically that the “complete” wave is well ap¬ 
proximated by merely the function i. With this result 
in mind, we are going, in the following, to derive exact 
analytic expressions for and for the (time domain) 
Green’s function in the complete SBH model. We shall 
then use these results to calculate the late time form of 
'I'l at null infinity in this model. 

The above iterative expansion appears to be an effec¬ 
tive calculation scheme for all modes I of the scalar radia¬ 
tion, except for the monopole mode 1 = 0. This is unlike 
the scheme used for the shell model in paper I, which 
held equally-well for all modes I with no exception. The 
reason for this difference between the two models in the 
monopole case will be discussed later. The calculation to 
follow shall regard only the modes with I > 0. 


IV. DERIVATION OF 

We first obtain an explicit expression for 'I'o) the first 
element of the iterative expansion. Since the only dis¬ 
continuity in the potential function Vq (at r* * = rg) is 
bounded in magnitude, we learn by the wave equation 
(0) that 4'g and its first order derivatives should be con¬ 
tinuous anywhere. In the sequel we explicitly use this fact 
in constructing an expression for ^'g. 

We shall consider separately three distinct regions of 
the domain S, as indicated in figure |^. Regions I (itg < 
u < —2rg) and II {u > —2rg,r* > rg) cover the part of S 
outside the surface r* = rg, while region III (r* < rg) is 
the portion inside this surface. 


i-i- 



FIG. 3. Calculation of th® complete Schwarzschild geom¬ 

etry. The solution in each of the regions labeled I, II, and III is 
discusses separately in the text. 


A solution to Eq. (0) for V = 0 has the general form^ 


^0/ = I] < 

n—0 


9oi [u) 

{v — uY~'^ 


(13a) 


'f'o// = 


A. 


(13b) 


'Soiii = F{u) + H{v), (13c) 


in which the labels I,II,III denote the region to which each 
specific solution corresponds, and where goi{u), goii{u), 
F{u) and H(v) are (yet) arbitrary functions. In the above 
equations the coefficients are given by 


, (2^-n)! 

" n!(;-n)!’ 


(14) 


and the parenthetical indices indicate the number of 
times the functions are differentiated. 

Causality implies that in region I, the solution 'I'g/)^) 
cannot be sensitive to the form of Vg at r, < rg. Thus 
it must be identical to the solution derived in the shell 
model, with the function goi{u) explicitly related to the 
initial data function by 

= (^) {u-uy-^T{u')du' (15) 


(see paper I for details). 

No w, with the initial condition '^oiii{v = 0) = C|^, Eq. 
( |l3c] ) implies that is a function of v only. We can 
then use the continuity of 'I'g_u at r* = rg to derive a 
closed differential equation for the function goji{u): 


This equation may be put into the form 
i+i 

^i?irg"[5ozz(«)](")=0, 

n—0 

in which the coefficient are given by 



(16) 


(17) 


^ The most general solution at r > ro involves also an arbi¬ 
trary function h{v). However, for our choice of initial setup 
(an outgoing initial pulse), the solution 4*0 can be expressed 
in terms of a function g{u) solely. This issue is discusses in 
appendix A of paper I. 

* Since we shall be interested mostly in the case where the 
initial pulse is specified at large distance, we assume here that 
its support is conhned to the exterior of the sphere r* = rg. 
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, 2-{2l-n)\ 

nl{l-n+l)l' 


(18) 


Thus goii (u) is a solution of a constant-coefhcients linear 
equation of order I + 1. It therefore admits the form 


die off exponentially with respect to their argument for 
all 1 < i < 1 + l|] We still have to derive an expression 
for 'I'o in region III, that is at r* < tq. By the continuity 
of 'I'o at r* = ro we have 'l'g^^(u) = 'l>g^(u = u —2ro). It 
follows that 


i+i 

9oii{u) = ^ Ci exp(-KiM/ro), (19) 

2=1 


/+1 

'^ 01 11 = ^ aiEi{v - 2ro), (26) 

i=l 


where Ci are constants, and the I + 1 complex numbers 
Ki are the roots of the algebraic equation 

i+i 

^i?i(-«)" = 0. (20) 

n=0 

The only properties of the numbers Ki important for our 
discussion are that (i) these numbers are all distinct (for 
any given value of 1), and that (ii) we have Re{Ki) > 0 
for all values of I and i. Hence go (and also dto itself) 
falls off exponentially at late retarded time u. 

The coefficients Ci are determined by imposing conti¬ 
nuity on ^'o at M = —2ro, namely by requiring (goii)^^^ = 
{goi)^^'^ at u = — 2ro, for all 0 < j < 1 . This leads to a 
set of Z -I- 1 algebraic equations for the I + 1 coefficients 
Ci , having the form 


i-i-i 

^ji^i = fo9oi 


i=l 


u— — 2rQ 


( 21 ) 


for 0 < j < I, where Mji = (—Kiy exp(2Ki). In a matrix 
form, we have 


det M = exp[2(Ki • -I- kz+i)] det K, (22) 
where K is the Vandermonde matrix 


/ 1 1 

— Ki —K2 


K = 


1 \ 


— K,\ —Kn 


V 




(23) 




which is always non-singular, provided only that the 
numbers Ki are all distinct (which is the case here). 
Therefore the set of equations (|^) has a unique solution 
for the coefficients Ci. 

W e finally obtain 'I'o// by substituting for goji in Eq. 
(13t). This yields an expression of the form 


I l-\-l / 

n=0 2=1 ^ ^ 


Ei{u), 


(24) 


in which are constant coefficients (being certain l- 
dependent functionals of the initial-data function r(it)), 
and where the functions 


Ei{x) = exp[-Kja;/ro] 


(25) 


where Ui = X)L=o Qmi atid where the functions Ei are 
those defined in Eq. (^). 

We conclude that 4>o “penetrates” the potential barrier 
only through a narrow null ray of typical width ~ 2ro 
adjacent to the initial ingoing ray u = 0. It has significant 
amplitude only in a “main” region ug < u ^ 0 and along 
that penetrating null ray. Elsewhere, 4/g is found to be 
exponentially small (in retarded time it at r* > rg or in 
advanced time i; at r, < rg) . This result (valid for all 
Z > I) is illustrated in figure |[ 
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FIG. 4. Domain of the function 4^0 Schwarzschild. Dark 
colored areas indicate regions where 4*^0 is not exponentially small. 
Note that 4^0 “penetrates” the potential barrier only through a 
narrow ingoing ray (of typical width ~ 2ro). 


V. CONSTRUCTION OF THE GREEN’S 
FUNCTION 

In this section we derive an analytic expression for the 
(retarded) Green’s function corresponding to the oper- 
ator dudy + Vo{r^.), with Vg{r) defined in Eq. (||). The 
(retarded) Green’s function G{u,v]u',v') is defined by 
the equation 


§ In Eq. (|^, as well as in all other expressions for the var¬ 
ious functions 4* to be appeared in this paper, it is to be 
understand that only the real (or alternatively—the imagi¬ 
nary) part is taken into account. The indication “Re” shall 
be omitted for brevity. 
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G,uv +VoG = d{u - u')S{v - v'), 


(29) 


(27) 

supplemented by the causality condition G{v < v') = 
G(u < u') = 0, where {u',v') are the null coordinates 
of a scalar “point” source (in the 1+1 dimensions rep¬ 
resentation), and {u,v) is where we evaluate the field 
this source induces. [It will become evident by construc¬ 
tion that this condition specifies a unique solution to Eq. 
(p^).] In view of the results obtained for 'I'o, we shall 
have to consider both ‘external’ (r(, > tq) and ‘internal’ 
{K < ro) sources. In what follows we treat each of these 
two cases separately. 


A. External sources 

We first consider a “point” source located at null 
coordinates {u',v') outside the surface r* = tq (thus 
v' — u' > 2ro). For this fixed source, we look for the 
Green’s function at any evaluation point (u,v). To that 
end we separate the future light cone of the point source 
into three regions, as indicated in figure Regions I 
and II correspond to evaluation points outside the sur¬ 
face r* = ro, while region III corresponds to internal 
evaluation points. 
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FIG. 5. Construction of the Green’s function for a scalar source 
sphere at r'^ > rg. The three regions I, II,III, defined with respect 
to that source, are treated separately in the text. 


We first observe that in Region I (that is at u < w' — 
2ro) the Green’s function cannot depend on the form of 
the potential at r* < ro (as implied by causality), and 
thus in this region is must be the same as in the shell 
model (outside the shell). Therefore, by Eqs. (33) and 
(35) of paper I we find that the Green’s function in region 
I reads 


G\u, V, u', v') = ^2 

n—0 


n 


{v — uy~'^ 


(28) 


in which the differentiation is with respect to u, are 
the coefficients given in Eq. (H), and 


ga{u-u\v') = 


l\ 


{v' — u){u — u') 


{v' — u') 


Now, in regions II and III Eq. ( |27[ ) is homogeneous, 
hence the solutions for the Green’s function in these two 
regions are of the form 




n—0 


{v — 


(30a) 


Qiii ^ (30b) 

where the functions g^^{u) and G^^^{v) are yet to be 
determined. 

By analogy with Eq. ( 0 ) we then have 
+1 

ffG = X! exp(-«;jMAo), (31) 

with Ki being the same numbers as in Eq. ([l^), and where 
the I + 1 coefficients Ci{u\ v') are to be determined such 
that the Green’s function is continuous along the ray u = 
v' — 2ro. This requirement leads to a set of ^ +1 equations 
for the coefficients Ci{u',v'), reading 

1+1 

= r^o[gGiu)y^^ 

(for 0 < j < /), where Mji = (—exp[—— 
2ro)/ro]. The solution (which always exists) is 

i 

Ci = exp[Ki{v' - 2ro)/ro] ^(■V^[g^(u)](^) 

u—v' — 2rQ 

j=0 



(33) 


with being the elements of the matrix reciprocal 

to the Vandermonde matrix (p^. Inse rting the explicit 
expression for g^ and using Eq. (30a), we can finally 
obtain for the Green’s function in region II, 


G 


II 


i 1+1 

=EEa 

n,j=0i—l 


{K-roy-^{roy+^- 


nji 


{<) 


iJ-n 


-Eiiu -v' + 2ro), 


(34) 


in which = (w' —u')/2, and where finji are certain con¬ 
stant coefficients (depending on I only). (Recall that the 
functions Ei die off exponentially with respect to their 
argument for all i.) 

To obtain the Green’s function in region III, we simply 
notice that G^^^{v) = G^^{u = v — 2ro), implied by the 
continuity of G at r* = rg. It follows that 


1 1+1 




j=0 i=l 


where /3ji = straightforward to ver¬ 

ify that with this result, we have G^^^{v = v') = I as 
necessary.) 
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B. Internal sources 

To obtain the Green’s function for a source point lo¬ 
cated at r' < To, we refer to figure |^, where again we 
indicate three regions, defined with respect to a given 
source at {u',v'). Again we discuss the construction of 
the Green’s function in each of these regions in separate. 
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FIG. 6. Construction of the Green’s function for a scalar source 
sphere at r* < rg. The three regions IV, V, VI are treated sepa¬ 
rately in the text. 


With Eq. (|^) , this constructs a set of / -I- 1 linear al¬ 
gebraic equations for the coefficients Ci. The solution 
reads 


Cj = exp[Kju7’’o], (39) 

where the numbers Ki are the same as for the external 
source. (Recall that the matrix K is always non-singular, 
hence this solution exist and is unique.) 

Using the results (^ and (^) we can finally obtain 


qVI 


l z-l-l 

n—0i—1 



l—n 

Ei{u - u'), 


(40) 


with the functions defined in Eq. (p5|), and where 'jni 
are certain constant coefficients (depending only on 1). 

To obtain G^, we simpfy notice that {v) = 

G^^{u = V — 2ro) (inferred by the continuity of the 
Green’s function), hence 


Z-(-l 

G'" =^7,£;7u-u'-2ro), (41) 

2=1 


where q* = EL=o 7ni- 


In region IV we have G^J(, = 6{v — v')8{u — u') by 
definition, which (by causality) leads to 

G^^ = 0{v — v')9{u — u'), (36) 

with 9 denoting the usual step function. 

In region V the Green’s function satisfies the homo¬ 
geneous equation G^^ = 0. With the continuity re¬ 
quirement G^(u = u' -|- 2ro) = 1, this means that 
qV _ qV Now, G^^ is given in terms of a func¬ 
tion gQ^{u), in a way analogous to G^^ in Eq. (30a). 
By the continuity of G_„ at = tq we must have 
= ^o) = 0, which is a linear differential equa¬ 
tion of order I + 1 for the function gQ^(u). The solution 
is (in analogy to Eq. (31), 


/ + 1 


5 g ^(«) = '^Ci{u', v') exp[-Kiu/ro], 


(37) 




with Ci being certain coefficients. 

To construct the coefficients Ci, we match the func¬ 
tion G^\ as inferred by Eq. (|3^), to its value on the ray 
u = u'. This value may be deduced independently by 
inserting the form G^^{u,v) = G^^{u,v)9{u — u') (im¬ 
plied by causality) into equation (p^, and observing that 
a solution must admit G,v = 0 along u = u'. This means 
that G is constant along this ray. By Eq. (^) (and re¬ 
quiring continuity) we then learn that this constant is 
unity. Requiring G^^{u = u') = 1 for all v then leads to 


[gjy]^^){u') = 0 
[g^Jf\u') = 1 


(0 < n < Z — 1) 


(38) 


C. Fixed external evaluation point 

Thus far we considered the Green’s function for a given 
sources at {u', v'), as a function of the evaluation coordi¬ 
nates (u, v). In practice, we shall be interested in calcu¬ 
lating the function at a given location (specifically— 
at null infinity, for M), what will involve integration 
over all possible sources. This requires knowledge of the 
form of the Green’s function at the evaluation location, as 
a function of the sources locations. To that end we only 
need to re-interpret our previous results: The expres¬ 
sions we have derived for the Green’s function shall be 
regarded as functions of the source coordinates {u',v'), 
with fixed evaluation coordinates {u,v). This reversed 
presentation of the results is illustrated in figure |^. Indi¬ 
cated in this figure are the regions of spacetime in which 
scalar sources influence the behavior of the scalar field 
at a fixed evaluation point (with null coordinates u, v) 
outside the surface r* = rg. Dark-colored areas in this 
figure indicate source regions where the Green’s function 
is not exponentially-small, as inferred by Eqs. d^), (|3^) 
and (p^). 
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FIG. 7. Construction of the Green’s function for a given exter¬ 
nal evaluation point at u, v. Regions where the Green’s function is 
not exponentially small are indicated by dark color. 

In region A {u' < u, u+2rQ < v' < v) Gis given by Eqs. 
( ^ ) and and is the same as in the shell model. In 
region B {v' < u + 2ro, r* > ro) G is given by Eq. (^. It 
vanishes here exponentially towards early advanced time 
v' , and it possesses significant amplitude only within a 
narrow ingoing null “band” (of typical width 2ro), adja¬ 
cent to region A. Then, in region C {u' < u, r'^ < rg) the 
Green’s function [given by Eq. (Effl] vanishes exponen¬ 
tially towards early retarded time^, and is of significant 
amplitude only within a narrow outgoing null “band” (of 
typical width 2ro) adjacent to the ray u' = u. 

Like in the shell model, we find that the main region 
of effective sources covers only the range u < v' < v 
of advanced times. In the shell model, however, the 
Green’s function vanishes identically outside this range 
(due to the divergent potential at the center of symme¬ 
try), whereas in the complete SBH model it dies off ex¬ 
ponentially (and also “penetrates” trough the finite po¬ 
tential barrier at u' = u). 

VI. CALCULATION OF vE'i AT NULL INFINITY 

We have shown that in the complete SBH model, d>o 
gives only an exponentially decaying contribution to the 
late time radiation. In this section we calculate the con¬ 
tribution of dti to this radiation at null infinity, and show 
that it is characterized by the same power-law tail of de¬ 
cay that was indicated in the shell model. Moreover, we 
show that even the amplitudes of the waves are the same 
on both models, provided that we choose [uqI ^ 2M ~ rg 
(the difference is of order rg/ug). 

Before we present the detailed calculation of dti, we 
first give some heuristic arguments concerning the ex¬ 
pected results. Figure ^ shows the region of spacetime 
in which scalar sources affect the behavior of the wave 
at null infinity, at a given retarded time u ^ M. Also 
shown, superposed, is the region where sources due to 
dtg exist. Outside the overlapping of these two areas, the 


Green’s function, or dtg, or both, are exponentially small. 
We expect (and later show analytically) that sources out¬ 
side the overlapping area shall give only an exponentially 
decaying contribution to dti at null infinity as u —> oo. 
We may thus focus only on the two overlapping regions 
shown in the figure. One of these regions lies inside the 
surface r* = rg (see the figure). It is of “dimensions” 
rg X rg, and is located near r* = (—u/2) <C M. In this lo¬ 
cation the potential function V (r* ) is exponentially small 
(see figure |^, and thus the contribution from this area 
should be exponentially small as well. We are left with 
the contribution of sources at the “main” region [namely 
Mg < m' < 0, It < r' < r], in which both the Green’s 
function and tkg are having the same form as in the shell 
model, except for in a narrow band (of width ~ 2rg) at 
the edge of this region. This suggests that for rg small 
enough, the calculation of should yield a result very 
close to that obtained in the shell model. 



FIG. 8. calculation of > 1/1 in the complete Schwarzschild model. 
The region filled with horizontal lines is where the amplitude of 
^0 is significant, while the region marked with vertical lines is 
where the Green’s function is significant, in accordance with the 
discussion given in the text. Potentially, significant contribution 
to is expected to arise only from areas where the two indicated 
regions overlap. 

To confirm the above heuristic indications, we shall 
now calculate at null infinity. In terms of the Green’s 
function derived above, we formally have 

nU nV 

d-5”(u) = - / dul dv' G°^{u]u',v')5V{u', v')^q{u', v'), 

J uq Jo 

(42) 

where 41^ and G°° stand for the value of 411 and G at 
null infinity (that is, for v oo). For the purpose of 
calculation, we separate the domain of integration into 
three regions. 
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poo pu pu-\-2rQ pv'— 2rQ pu-\-2rQ pu 

/ dv' / du' + / dv' / du' + / dz;' / du' (43) 

n+2ro *^0 *^140 0 Jv'—2rQ 

to be labeled A, B, and C, respectively, as indicated 
in figure In what follows we consider separately the 
contribution from each of these three regions to 4'“. We 
show that the contribution of region A is the dominant 
one, and that this contribution is identical (in its pattern 
of late time decay, and also, to a certain order of accuracy, 
in its amplitude) to that obtained in the shell model. 


Contribution from region C 


In terms of the new integration variables = {v' — 
u')/2 and t' = {y' + u')/2, and using Eqs. (^) and (^ ) 
(with r* ^ oo), the contribution to from sources in 
region C takes the form 


i-l-l 

pro 

pu-\-r^ 


/ dr'. 

/ dt'Ei{u — t' + r',) 


J —u/2 J 

— r'^ 

X Vir',)E,it'+r', 

- 2 ^ 0 ). 

(44) 


(recall that in region C, where < rp, we have 5V = V 
by definition.) 

Since V(r'^) oc exp[(r(/2M)] for ^ —oo, one finds 
that the integrand in the last expression dies off expo¬ 
nentially in retarded time u anywhere inside the domain 
of integration. It is easy to verify that the integral it¬ 
self would be exponentially small for large u/M. For 
example, for any fixed retarded time u ^ M there exist 
positive constants ci, C 2 , C 3 , C 4 and k, such that the fol¬ 
lowing upper bound is applicable to the above integral 
(in absolute value): 


Z +1 


r*—n/4 








u/2 




+ C 2 dr'^ / dt'\Ei{u - t'+ r'^)\ 
J — uj A. J — r'^ 

pTQ 


+ C3 y dr'^ y^^^, dt' \Ej{t' -I- r(, - 2ro)| | 

< C 4 exp[—«:(u/ro)]. (45) 

We conclude that internal sources of 'I'o (namely 
sources at r* < tq) give at most an exponentially decay¬ 
ing contribution to the late time radiation at null infinity. 


Contribution from region B 


By Eqs. (13b) and (|^) (with r* ^ 00 ), the contribu¬ 
tion from region B to 'I'l at null infinity reads 


i+i i 

»fB = E , 

2—1 0 ^ 

r'- 2 ^o(r(-ro)'-^r^o 


u+2ro 

njii dv'Ei{u - v'+ 2ro) 


du 


iO 


21 —n 


SV{r:)glr\u'), (46) 


in which are certain constant coefficients, and the 
function go(u') stands for the expressions derived in the 
previous chapter for gl{u') and {u') [Eqs. (H and 
(|T^), respectively]. If we now integrate this expression 
by parts with respect to v', we find that to the leading 
order in M/u and in vq/u, 


i i 

n—Oj—0 


^TB = y:y'jnjro 


{u — u'Y ^rl 


X 6V{u — u + 2ro) g^yu) 


/«o - w'+ 2ro)^' ” 


(47) 


with Pnj being some other coefficients. Now, integrate in 
parts each of the terms n successive times with respect 
to u'. The resulting surface terms would all be negligible 
at large u/rp, since go{u') dies off exponentially at large 
retarded time u [see Eq. (p^j. In addition, these surface 
terms are strictly compact from below. We are left with 


i i 


^^(-l)"/3„,ro / 

n =0 4 = 

0 

d” 

(u — u'y~h 

du'^ 

{u- u' + 2ro) 


21 


—SV {u — u' + 2rH.) 


(48) 


to the leading order in M/u and in r^ju. 

To continue, we shall have to write 5V in terms of the 
null coordinates. This cannot be done explicitly, since 
the function r(r*) is implicit. Rather, we shall use the 
large r expansion 

dV(r^ > ro) = Mr~^[a -\- 51n(r*/2M)] 

+ &{M‘^rZ^[\n{r,/2M)f) , (49) 

where a and b are constant coefficients, depending only 
on I and M. In paper I we argued that (in the framework 
of the shell model) it is merely the asymptotic form of 
the background potential which affects 4'“ at u M. 
This has been also tested numerically (see figure 11 in 
paper I). We now proceed by assuming that the same is 
true in the complete SBH model as well. 

With SV taken to the leading order in M/r^^ Eq. ( |4^ ) 
takes the form 


I )_i j^U / !\ 




j=0 m=0 •^“0 


{u-u' + 2ro)2'-™+3 


^jm 4“ bjjnln 


u — u' -\- 2 ro 


(50) 


to the leading order in M/u and in r^/u^ where Ujra and 
bjm are certain constant coefficients. We observe that, 
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since go{u') ~ exp(—u'/ro) at large u', the upper part of 
the integration over u' (say, between u' = '/uM and u' = 
u) gives a contribution which dies off exponentially at 
large u. In our approximation we can thus concentrate on 
the contribution coming from early retarded times (say, 
uq < u' < VuM)- At any large enough retarded time 
u there exist positive constants C 5 and cq such that this 
early contribution (in absolute value) is bounded from 
above by 


where akj and b^j are certain constant coefficients that 
depend on I and M, but not on tq. Integrating over v', 
the r. h. s of the last equation becomes 


l l-k k 


^EEE 


, (u - uy-^-^ro f-^ 

du T-: — TTTT —^ -TT X 


k—O 3—0 m —0 

dkjm H“ ^kjm ln 


(u - u' + 2ro)2'-J-™+2 
u — u' + 2ro 


2M 


go{u'). (55) 


i i-j 

|^'?bI < CsMro EE 

j—O m—O 


— uo]{u — uoY ^ "“tq 


[u — 


ffi l^jml la 


u — Uq + 2ro 


2M 




(51) 


In what follows it will become apparent that this contri¬ 
bution to the late time radiation at null infinity dies off 
more rapidly than the radiation due to scattering in the 
“main” region of sources (region A) , which will be shown 
to be characterized by a decay tail. Therefore, the 

contribution from region B is negligible at u ^ M. 


Contribution from region A (the “main” region) 

The remaining contribution to calculate is that coming 
from region A, reading 


with Okjm and b^jm being yet another constant coeffi¬ 
cients, independent of rg. 

Now, since go{u') falls off exponentially at large u', we 
may cut off the integration at, say, u' = without 

affecting the integral to the leading order in M/u. Doing 
so, we observe that the leading order contribution to this 
integral comes only from terms corresponding to m = 
k. (Note the way the dependence in the parameter rg 
cancels in the leading order). 

Defining 

pVuM rOO 

/ goi{u)du' ~ / goi{u)du' = /g, ( 56 ) 

J UQ ^ —00 

(were the first equality holds to the leading order in M/u, 
as go is compact from below and dies off exponentially 
at large u) we find that to the leading order in M/u, in 
rg/u and in ug/u, 

= M/gu-'-2 [fci-ffc 2 In u/M]. (57) 



To evaluate this expression we first write d /g(u' , v') in 
terms of a function go{u') [as in Eqs. (|l^a),( |l3^ ], then 
integrate by parts each of the resulting terms in the r. h. s 
n successive times with respect to u'. Neglecting surface 
terms, which are all exponentially small at late time since 
go ~ exp[— u/rg] at large u, we obtain 


^ pU nOC 

^?a(«) =-E-l)”^" / / dv’goiu') 

„_g Juo Ju+2ro 


9u'” 


G°^{u-,u',v')5V{u',v') 


{v' — 


(53) 


With the explicit form of the Green’s function [Eqs. 
( p^ and (29)], and with SV taken to the leading order 
in M/r^:, the last equation takes the form 


l l-k 


I -r» Un *• 


dv 


, {u-uy-^-yv' -u) 


k=0j=0''^0 ■'u 


+ 2 ro y - u'y-:>+^ 


X Clkj H” ^kj In 


2M 


go{u'). (54) 


ki and ^2 are constant coefficients that do not depend on 
rg. The only remaining reference to the value of rg lies 
within the integral Jg. 

Since the values of the coefficients ki and ^2 are inde¬ 
pendent of rg, then in order to obtain these values one 
may use Eq. (m) with whatever value of this parameter 
(requiring only that rg ^ u). Now, if we take rg = 0, 
then Eq. (p^ ) becomes completely analogous to the ex¬ 
pression for 4'/° in the shell model (Eq. (37) in paper 
I). Comparing these two expressions, we find that the 
Green’s function G and the potential 5V appearing in 
both integrands are exactly the same. The two expres¬ 
sions differ only in the form of the function 'kg, whi ch o n 
both cases is expressed in a similar way [as in Eq. (13a)] 
in terms of two different functions go{u)- However, the 
explicit form of the function go{u) (as related to the ini¬ 
tial data) has no effect whatsoever on the value of the 
coefficients ki and k 2 and so these coefficients must be 
the same as in the shell model. Therefore, comparing 
Eq. d^ ) with Eq. (52) of paper I, we learn that ^2 = 0, 
and that, to the leading order in M/u, in ug/u, and in 
rg/u we have 


4'/“(u » M) = 2{-iy+\l + ly.MIo 


- 1-2 


(58) 


[for the compact initial data setup]. 

We conclude that the wave has the same late time 
behavior in the complete SBH model as it has in the 
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shell model, namely it dies off as provided that 

the initial pulse is compact. Numerical calculation of 
in the complete SBH model agree with this result, as 
demonstrated in hgure 



FIG. 9. Late time tails of at null infinity. Presented on 
a log-log scale are numerical results obtained for in the com¬ 
plete SBH model at u = 10^M (approximating null-infinity), for the 
/ = 0,1, 2, 3 modes. (For I = 0 we used the definition of given at 
the end of section VI). Compact initial data for the numerical prop¬ 
agation has been specified between u = — 40M and u = — 50M, and 
the parameter tq has been set to 3M. Also shown, for reference, are 
dotted lines proportional to u~^~^. The results demonstrate the 
u~^~^ late time decay rate predicted by the analytic calculation. 

We may similarly obtain the u~^~^ decay characteriz- 
ing the static initial setup, by comparing Eq. (^) with 
Eq. (54) of paper I. If a static scalar field is present out¬ 
side the central object up to some moment before the 
event horizon forms (no static solution exist which is well 
behaved both at the event horizon and at inhnity), we 
shall have, to the leading order in M/u, in uq/u, and in 
ro/u, 

vI/-(u»M) = (-2)'+i|^Ma.u-'-\ (59) 

where /i represents the amplitude of the initial static held. 

It is also instructive to compare the amplitudes of the 
wave at some hxed u ^ M value on both models 
(the shell and the complete SBH), given the same initial 
data. As implied by the above discussion, the relative 
amplitude is simply given by the ratio of the integrals 
Jo = go(u')du' associated with both models. Relying 
on the explicit expressions derived for the functions go in 
this paper and in paper I, it can be easily shown that the 
two integrals Jq, associated with the two models, differ 
merely by an amount of order ^ tq/uo- Thus, concen¬ 
trating on the case |uo| 2M ~ tq (in the context of 
which our analysis proves effective—see the discussion in 
the following section), one observes no difference between 
the late time behavior of the wave 'I'l at null infinity on 
both models. This result is accurate to the leading order 
in M/u, in uq/u and in to/uq- In particular, we may 


conclude that at late time the wave has no refer¬ 
ence (in our approximation) to the value of neither the 
radius of the shell R (in the shell model) nor the param¬ 
eter tq (used in the complete SBH model). This result is 
consistent with the assumption that details of spacetime 
structure at small r values are not manifested in the form 
of the late time radiation. 


The monopole case (I — 0) 

It was pointed out while introducing the iterative ex¬ 
pansion, that a scheme based on that expansion fails to 
handle correctly the case 1 = 0. In the monopole case 
it is straightforward to find that a Green’s function de¬ 
fined as in sec. ^ is constant (a unity) throughout the 
whole range of evolution. It then follows that the wave 
'I'l is constant at late time, resulting in the divergence of 
higher terms of the expansion (d> 2 , '£' 3 ,...). Looking for 
the cause of this failure, we notice that in the cases I > 0 
it is the centrifugal potential barrier that “cuts off” the 
Green’s function and confines it (for late retarded time 
evaluation points, u M) mainly to late retarded times 
{u < u' < v). In the shell model (see paper I) it was 
the presence of the center of symmetry which effectively 
acted as a potential barrier for the Green’s function even 
in the monopole case, where no centrifugal potential ex¬ 
ists. This is why the iterative expansion applied in the 
framework of the shell model proved to be equally effec¬ 
tive for all modes of the radiation. 

To analyze the case Z = 0 in the complete SBH model, 
one is thus led to try a different iterative expansion, de¬ 
fined such that the Green’s function is subject to an ap¬ 
propriate potential barrier, as for the modes with I > 0. 
One technically simple possibility it to take 

Vt° = M-^5{r,). (60) 

We then define the iterative expansion as in Eqs. 

with the ‘new’ potential V/~°. With this def¬ 
inition, we find that = r(u), and that the Green’s 
function (at an evaluation point u,v with (v — u)/2 = 
r* > 0 ) is given by 

G{u' <u,u<v' <v) = l 
G(u'< u'< u) = exp[(u'— u)/M] . (61) 

g\v' <u'<u) = exp[(u' — u)/M] 

(The three regions indicated in this equation are those 
labeled ’A’,’B’, and ’C’, respectively, in Fig. (|^), when 
setting ro = 0). A simple calculation [based on Eq. (|4^)] 
then shows that at null infinity is given by Eqs. 

( ^ ) and d^ ) (with 1 = 0). These equations are therefore 
valid for all modes 1. 
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VII. HIGHER TERMS OF THE ITERATIVE 
EXPANSION 

The Green’s function technique applied thus far pro¬ 
vides a formal way to calculate each of the terms n one 
at a time, in an inductive manner. We shall not, how¬ 
ever, calculate further terms of the expansion, but rather 
we will refer to the results of the analysis in paper I. 
Regarding the shell model, strong indications were given 
that 

• The dominant contribution to ^ M), for all 

A > 1, is only due to sources of at large 

distances. 

• All terms of the iterative expansion (excluding 

'I'o) seem to share the same late time pattern of de¬ 
cay at null infinity, namely a inverse power- 

law tail (for compact initial pulse), or a u~^~^ tail 
(for static initial held). 

• If the initial pulse is conhned to large distance, 
|uo| 2> M, then the iterative sum converges at null 
inhnity rather efficiently to the ‘complete’ scalar 
wave 'h at late retarded time u. 

• Moreover, in this case (|uo| ^ M), the scalar wave 
d/ is well approximated by merely ^*1 (with correc¬ 
tions smaller by order M/uq). 

Now, it is reasonable to assume that all four of the 
above results are also valid in the complete SBH model. 
For, in both models (the SBH and the shell models) 
spacetime structure at large distance is the same, and, 
as argued in paper I, it is the large distance region whose 
structure is relevant in determining the late time form 
of the waves at null-inhnity. (This conclusion has been 
demonstrated in the framework of the shell model by ex¬ 
plicit analytic calculation of 'I'l and ^*2 at null-inhnity. 
Physical explanation was given in the concluding section 
of paper I. For the complete SBH model, a demonstration 
is provided by the explicit analytic calculation of in 
the preceding section). Actually, we need only to assume 
that the hrst of the four results indicated above holds 
in the complete SBH model. Then, the same reasoning 
used in section VH of paper I in deriving Eq. (72), leads 
us immediately to realize that ~ in the SBH 

model as well. Also, a completely analogous analysis to 
that applied in section VHI of paper I shows that the 
third and forth of the above results hold in the complete 
SBH model as well. 

Numerical analysis of the complete SBH model hrmly 
supports the above arguments, showing that all four re¬ 
sults are indeed valid in the complete SBH model. In 
what follows we present some examples of these numeri¬ 
cal experiments. 

Figure |I^ presents the ratio calculated nu¬ 

merically for the monopole and the dipole modes, for 
various values of the parameter uq. (The “complete” 


wave d/ has been obtained by a direct numerical solution 
of Eq. (H)). The results demonstrate that (like in the 
shell model) as \uo\/M is set larger, 'hi becomes a better 
approximation to the “complete” wave 'k at null infinity 
at late time. 

In hgure ^ it is demonstrated numerically (for I = 
0 ,1,2) that the iterative series applied in this paper seems 
to converge rather efficiently for a large \uo\/M value. 



FIG. 10. ' 1^1 approximates the “complete” wave at null in¬ 
finity at late time, provided that the parameter \uo\/M is chosen 
large. This is demonstrated numerically for the monopole {I = 0) 
and the dipole {I = 1) radiation by comparing the ratios at 

null infinity (approximated hy v = 10^ M) for the various values 
(a) UQ = —40M, (b) UQ = —200M, (c) uq = —lOOOAf, and (d) 
UQ = —5000M. The parameter ro is set to 3M. 
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u/M 

FIG. 11. Numerical indications for convergence of the iterative 
scheme at null infinity. Presented (on a linear scale) are the ratios 
(a) 'I'l/'I', (b) and (c) sample values 

Z = 0,1, 2. The other parameters are set to rp = 3, uq = —200M 
and V = 10®M (approximating null infinity). The results suggest 
a rather efficient convergence of the iterative expansion for large 
\uq\/M values at large retarded time u. 


VIII. TAILS AT CONSTANT RADIUS: THE LATE 
TIME EXPANSION 

So far we were discussing an analytic technique en¬ 
abling the calculation of the late time behavior of the 
scalar field at null infinity. In this section we apply a dif¬ 
ferent, local, analysis to study the late time behavior of 
the wave along any r = const world-line outside the black 
hole, and along the event horizon. Using this method we 
will be able to derive a simple analytic expression for the 
field, consistent with an inverse power-law decay, and 
accurate to the leading order in M/t (or in M/v along 
the horizon). However, this expression will inhere two 
undetermined parameters (one for the power-law index 
and the other for the amplitude). We shall deduce these 
parameters by matching our late time solution at null 
infinity to the form derived in the previous sections us¬ 
ing the iterative scheme. In that respect, the iterative 
scheme shall prove to be an essential key for the con¬ 
struction of a complete late time description of the wave 
behavior anywhere outside the black hole. The purpose 
of this section is to introduce and apply the local method. 

We make the assumption that at late time, the Klein- 
Gordon scalar wave [priorly SOTarated in terms of the 
spherical harmonics as in Eq. (g)] admits the expansion 

OO 

0 * (r, v) (62) 

fc =0 

in which the number ko and the set of functions Fk (r) are 
yet to be determined. Substituting this expression in the 
Klein-Gordon equation (^ and collecting terms of com¬ 
mon u-power, the partial equation is thereby converted 
to an infinite set of ordinary coupled equations for the 
unknown functions Fk(r), 

^2 ^ p ^ 

2(fco + A:-l)(r2F'_i+rFfc_i) (63) 

(for fc > 0), where a prime denotes differentiation with 
respect to r, and where we have set F_i = 0. This set 
of equations exhibit only a ‘weak’ coupling, in the sense 
that each of the functions Fk depends only on its preced¬ 
ing function Fk-i, with Fq obeying a closed homogeneous 
equation. This hierarchy allows one to treat each of these 
equations one by one, in an iterative way. In this pro¬ 
cedure, each of the functions Fk (with fc > 1) satisfies a 
closed second-order inhomogeneous equation. 

We proceed as follows: First, we show that there exists 
a solution (j)^ of the form (|6^, which is regular anywhere 
outside the black hole, in particular at the event horizon 
and at null infinity. Then, with the aid of our previous 
results at null infinity, we deduce the late time behavior 
of the scalar waves as detected by a static observer at any 
constant radius. In particular, we obtain the late time 
form of the waves along the event horizon. 
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where ak and bk are arbitrary parameters, and 


We define a new dimensionless radial coordinate 


P = 


r — M 
M 


(64) 


which varies monotonically from the event horizon (p = 
1) up to space-like infinity (p = oo). In terms of the new 
variable, Eq. (|^) takes the form 

(p" - 1)F'' + 2pF' - l{l + l)Fu = Dfc(p)[Ffe_i], (65) 


in which a prime now denotes differentiation with respect 
to p, and where D[p) is the differential operator 


Dk{p) — 2M{kQ + k — l)(p -l- l)[(p + l)9p -l- 1 ]. ( 66 ) 

(Note that Dk{p) depends on fcojbut is independent of 
1). For fc = 0, the r. h. s of Eqs. (|^ ) vanishes. 

We would like to construct solutions Fk to Eqs. (^) 
such that (j) would be regular both at the event horizon 
and at infinity. To allow (j) to be regular at the horizon 
(where the v coordinate takes finite values), we require 
all functions Fk (p) to be regular at p = l[]. We will now 
show that such regular solutions Fk{p) do exist^ 

We shall construct the functions Fk in an iterative 
way, starting with Fq. For k = 0, Eq. ( [^ is ho¬ 
mogeneous, and its general solution is given by Fg = 
aoPiip) + boQi{p), where ag and bg are arbitrary param¬ 
eters, P/(p) is the Legendre polynomial of order I, and 
Qi (p) is the Legendre function of the second kind, of or¬ 
der 1. The polynomials Pi{p) are, of course, finite at the 
event horizon (p = 1 ) and divergent (as p*) at p —*■ oo. 
Conversely, The functions Qi (p) diverge at the event hori¬ 
zon and vanish (as p~^~^) at p —> oo. Regularity of Fg 
at the event horizon therefore requires that 6 = 0 , hence 
we obtain 


F'(p)=aoFz(p). (67) 

Now consider Eq. (|^) for a general function Fk (with 
k > 1). The general solutions to the inhomogeneous 
equations read 


FUp) = akPiip) + bkQiip) 

fP Qi{p')Dk{p')[Fk-,{p')] 


Piip) 


- Qi{p) 


(p'2 - l)W(p') 

Pi{p')Dk{p')[Fk-i{p')] 

(p'2 - l)W(p') 


dp' 

dp', ( 68 ) 


** Had we apply a 1/t expansion of 4> instead of the 1/v 
expansion (^), all functions Fk had to diverge at the horizon 
(where t = oo) to assure regularity of the scalar wave there. 
For that reason, the 1/v expansion seems more plausible from 
the technical point of view. 

^Tn Mathematical terminology, Eq. (|65|) possesses a “regular 
sing ular point” at p = 1 (the event horizon)—see, for example, 
[|l5[ . For this case, standard mathematical theory tells that 
series solution for Fq and for Fi can always be constructed 
near p = 1. However, for k > 2 the source term in Eq. 
( psj ) appears to involve logarithmic functions, for which case 
standard theory gives no clear rules. 


W = P/Qi - PiQ'i = (p2 - 1)-1 (69) 

is the Wronskian. Using the relation 

Qi{p) = Pi{p) ^^^dp' (70) 

and integrating Eq. (^) by parts, we can then obtain 

Pki.P) = akPiip) + bkQiip)+ 

Pi{p) j^dp' ^j^J) dp"Pi{p")Dk{p")[Fk-iip")]. (71) 

We now show by mathematical induction that with 
bk = 0 (for all k), the functions Fk{p) are all analytic 
at the event horizon. The first function, Fg, is analytic 
at p = 1 by Eq. (|^. Now, following the inductive pro¬ 
cedure, assume that Fk-i is analytic at p = 1 for some 
fc > 1. Then, Dk{p")[Fk-i{p")] is analytic at p" = 1, 
hence the integrand of the p" integration in Eq. ( 0 ) is 
analytic at that point. We thus find that the integral over 
p" can be written in the form (p' — l)/(p'), where f{p) is 
some function which is analytic at the horizon (this can 
be shown by expanding the integrand in a Taylor series 
near p" = 1 , where it is analytic). Since the polynomials 
Pi have no real zeros is the range p > 1 , and W 

diverges as (p' — 1 )“^ at the horizon, we conclude that 
the whole integrand of the p' integration is analytic at 
p' = 1 , and therefore that the integral over p' must be 
analytic as well. Hence the solutions Fk, defined in an 
inductive way by Eq. with 6 fe = 0 for all A: > 1, are 
all analytic at the horizon. 

By this we have shown that the wave equation admits 
solutions (p of the form (|^), which are analytic at the 
event horizon. The most general of these solutions con¬ 
tains an infinite number of free parameters, one for each 
power of 1 /v 

We do not know yet the value of the power index fco, 
appearing in the expansion (p^. To obtain this value we 
shall now evaluate p at null infinity. By this mean we 
will be able to (i) show that the form of (p at null infinity 
is consistent with the results of our iterative analysis (in 
particular, that p is regular there); and (ii) deduce the 
value of fcg by comparing the results arising from the two 
independent schemes. 

We start by showing, using mathematical induction, 
that the functions Fk all have the asymptotic form 

Ffc(p ^ oo) ~ aoCfep'+'', (72) 

in which Ck are certain constant coefficients, other then 
zero and yet to be determined. Here and henceforth, the 
form /(x) ^ ex'" (where c is some constant) means that 
lim 2 ;^oo[/(x)/x"] = c. It appears most convenient to 
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(73) 


{2l)\{2l + l)U 

(Z!) 2 (fco-l)!' 


(81) 


prove Eq. ([^ by first showing that[^ 



oo) ^ ao{l + k)ckp’'~^^ 


Then, Eq. ( [^ i s implied. 

The form ([73]) obviously applies for Fq [given in Eq. 
(|^)], with 


(2; + l)!! 

d( 2 ; + 1 )’ 


(74) 


which is just the coefficient of p* in the polynomial Pi{p)- 
Following the inductive procedure, we now assume that 
Eq. (^^ applies for some fc > 0, and show that this leads 
to dFk+i/dp ~ (/ + fc + l)aoCfc+ip*+*. Our assumption 
necessarily implies that Fk ^ a^Ckp’’^^- Hence, by Eq. 
( 1 ^ ) we have 

Dk+iFk ~ 2MaoCkiko + k){l + k + l)p^+^+\ (75) 


Consequently, for the integration over p” in Eq. & we 
obtain the asymptotic form 


2MaoCkCo 


(fco T k){l + /c + 1) 


{2l + k + 2) 
ft then follows from Eq. (0) that 


ip‘ 


/ \ ‘2l-\-k-\-'2 


Fi 


fc+l ^ OoCfc+lP 


i+fc+1 


with 


Cfc+i — 2M Ck 


(fco + k){l + fc + 1) 

(2; + fc + 2)(A: + 1)' 


(76) 


(77) 


(78) 


Finally, differentiating Eq. we get dFk+i/dp {I + 
k + l)aoCfc+ip*+^, which establishes the inductive proof 
ofEq. (||). 

We have thereby shown that the functions Fk all admit 
the asymptotic form (|7^), with the coefficients Ck given 
by the recursive formula (|^), supplemented by Eq. (lA). 
We obtain, in conclusion, 


Fk^aoai{2fCkM-Wi+^ (79) 


for r* ^ oo, where we have explicitly used the fact that 
p ~ r/M ~ r^,/M, and where 


{ko + k - iy.{i + k)\ 

{21 + k + i)\{ky. 


(80) 


and 


(The coefficients Ck and ai are not to be confused with 
the coefficients appearing in sec. (p^)). 

Eq. ([^ describes the form of the functions Fk to the 
leading order in r/M, which is sufficient for our purpose: 
matching (/ at null infinity. We comment, however, that 
a full series expression for the functions Fk at large r can 
be obtained as well. It has the form 

k 

Fk{r) = Y.^rlMy-^Hk,{M/r)Mr/M)Y, (82) 

j=o 

where Hkj{M/r) are Taylor series. This form can be 
verified by substituting it into Eq. d^), then constructing 
explicit recursive formulae for the coefficients of each of 
the various series Hkj. Note that to the leading order 
in r/M no logarithmic terms are involved, and the form 
(^) is recovered. 

To obtain </> at null infinity, we insert Eq. (^ into the 
expansion (^). We get (to the leading order in r/M), 

= aoai{2M)-‘v^-^° II (l - . (83) 

k=0 

To evaluate the power series, we write it in terms of a 
generating function. 






^- 21-1 


^ko-2l-2 

dqko-2l-2 


qko-1 


yL 

dq^ 



, (84) 


which is valid for \q\ < 1. In this expression the deriva¬ 
tives might be of negative orders, in which case inte¬ 
grations are implied. If we now make the substitution 
q = (l — ^), keeping just the leading order in u/v, we 
find that 


J2Ck{i 


k=0 


U\k _ / 1 \ 

v) dq^<>~’-~^ \1 —< 7 / 


{ko-l 



ko —l — l 


(85) 


Since A! = rcj) = vcj) to the leading order in u/v, we finally 
obtain 


^-“(u) = aoaiiko - I - 2y.{2M)-^u-^^°-‘-^\ (86) 

where 'k°°(u) stands for the wave T evaluated at null 
infinity. 

The value of both the power index ko and the yet-free 
parameter ag can now be specified by comparing the last 
result to the results arising from our iterative scheme, 
Eqs. (1^ ) and (^^. This comparison yields (assuming 
that juoj ^ M), 


kn = 


«This is not valid when I = 0 and fe = 0, for which case 
dFk/dp = 0. However, the rest of the analysis does not 
change. 
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2^ -I- 3 no initial static field 
2^ -I- 2 initial static field 


(87) 
















and 


Oo 


(V.f(2l+2)\ , 2ilfy+^/ri 
(2Z)!(2i+l)!! V -'0 

(i!) (2i+l) / lAfV+lu 
2(20!(2Z+1)!! V 


no initial static field, 
initial static field 


( 88 ) 


where the integral Iq [defined in Eq. is directly re¬ 
lated to the initial data via Eq. (HI), and /r is the ampli¬ 
tude of the initial static field (when present). 

Provided with an exact expression for Eo(r) and with 
the value of fco, we are now in a position to write the 
form of the scalar field at any hnite value of r, at very 
late time. By Eq. ( |6^ we have, to the leading order in 
M/t and in M/uq, 


= (i»k*l); (89) 

and to the leading order in M/v and in M/uq, 

event horizon, (90) 

where the two values of powers correspond to the cases 
where an initial static field is or is not present, respec¬ 
tively. 

We should emphasize here that the above results [Eqs. 
( ^ ) and (0] apply towards time-like infinity at any 
value of r, establishing [together with the results at null 
infinity, Eqs. (p^, (|^)] a complete picture of the late 
time behavior outside the black hole. To the best of our 
knowledge, such a result has never been obtained pre¬ 
viously. (For example, in and Q, analytic ex¬ 

pressions were derived only for the asymptotic domains 
r* » M and r* <C M.) 


IX. CONCLUDING REMARKS 

In this paper, and in the paper preceding it, we have 
tested the applicability of a new analytic scheme for the 
calculation of the late time behavior of fields outside 
black holes. It was demonstrated, considering the simple 
model of scalar waves outside a SBH, that a simple ex¬ 
pansion of the held near time-like inhnity can be used in 
order to construct a late time solution consistent with a 
power-law decay anywhere outside the black hole. How¬ 
ever, the actual index of the power-law, as well as the 
amplitude coefhcient of the wave (as related to the ini¬ 
tial data) could not be determined merely by this local 
analysis. This information could be obtained only by a 
full integration of the 2-dimensional initial value prob¬ 
lem for the wave evolution, technically enabled by the 
introduction and application of the iterative procedure. 

Thus, by applying both the iterative scheme and the 
late time expansion, we were able to obtain an analytic 
expression for the scalar held in Schwarzschild, accurate 
to the leading order in M/t (or in M/u at null inhnity, 


or in M/v at the event horizon) and holding anywhere 
outside the black hole. The expression calculated is ex¬ 
plicitly related [via Eqs. (15), (^), and (0] to the form 
of arbitrary initial data specihed at large distance (our 
approximate solution has corrections of order M/uq). 

Of course, the main justihcation for the introduction of 
the new approach should rely on its applicability to more 
realistic models of wave evolution, for which no other an¬ 
alytic approaches have been proposed. In what follows 
we mention some possible applications of our calculation 
scheme, which include the analysis of scalar helds outside 
Kerr black holes, the analysis of gravitational perturba¬ 
tions in Kerr, and the extension of our analysis to the 
interior of black holes. 

The most interesting application of the new scheme 
concerns rotating black holes. As already mentioned in 
the introduction, this generalization is the prime motiva¬ 
tion behind the presentation of our approach, since real¬ 
istic stellar objects (and black holes) generically posses 
angular momentum. The generalization of our analysis 
to the case of scalar waves propagating in the exterior of a 
Kerr black hole shall be presented in a forthcoming paper. 
In brief, the basic idea behind this generalization is to ex¬ 
press the lack of spherical symmetry in Kerr spacetime in 
terms of interactions between the various modes of spher¬ 
ical harmonics. The resulting interaction terms coupling 
the field equations for the various modes (these terms 
are expected to be small, in a sense, at late time) are 
then to be treated using our iterative technique. Apply¬ 
ing an iterative decomposition basically similar to that 
used in the spherically-symmetric models, those interac¬ 
tion terms become source terms in the resulting hierarchy 
of wave equations. The mathematical treatment of these 
equations is then similar, in principle, to that applied in 
the spherically-symmetric cases. This provides the late 
time form of each of the modes at null infinity. Then, a 
generalization of the late time expansion method (based 
on the same interaction-between-modes approach) pro¬ 
vides the late time behavior of the field anywhere outside 
the Kerr black hole, in particular along its event horizon. 
The details of both parts of the analysis in Kerr shall be 
given in (see also ED- 

Obviously, the scalar model discussed so far is just a 
simplified analogue to the physical problem concerning 
the dynamics of gravitational perturbations. The plausi¬ 
bility of the scalar model stems from the remarkable re¬ 
semblance of the underlying mathematical formulation, 
between this model and realistic models of gravitational 
waves (as was already realized in [|^, for example, for the 
SBH case). Equations governing metric perturbations of 
the SBH were derived by Regge and Wheeler 0 (for 
axial perturbations) and by Zerilli ||l8|| (for polar pertur¬ 
bations). Both equations can be put in the same form 
as the scalar field equation (|^) , where this time the wave 
function represents certain linear combinations of entities 
characterizing the metric perturbation. In both Regge- 
Wheeler’s and Zerilli’s equations one also finds that the 
effective potential is similar in shape to that of the scalar 
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model (Eq.(6) and figure 1). This suggests that the prob¬ 
lem of gravitational waves propagation in the SBH ge¬ 
ometry may be treated using the same iterative scheme 
applied for the scalar model. 

We do not have similar separable equations for met¬ 
ric perturbations of Kerr black holes. Rather, a sec¬ 
ond approach, based on the Newman-Penrose tetrad for¬ 
malism, was used by Teukolsky to derive separa¬ 
ble wave equations governing perturbations of the Weyl 
scalars. [An exhaustive discussion of field equations for 
gravitational perturbations is given in |20|, pp. 174-182 
(for Schwarzschild) and pp. 430-443 (for Kerr)]. Sep¬ 
aration of Teukolsky’s equation is only possible in the 
frequency domain (namely by first separating the wave 
to its Fourier modes). To apply our iterative approach, 
however, the time dependence should rather be kept in 
the master Teukolsky equation. Instead, the master per¬ 
turbation equation is to be treated in the way outlined 
above, namely, by considering interactions between the 
various modes of spherical harmonics. By now we have 
first indications that this approach is indeed applicable 
to gravitational perturbations (in the tetrad formalism). 
We intend to study this subject more deeply in the near 
future. 

Finally, we mention the possibility of extending our 
analysis to internal perturbations of black holes. Re¬ 
cently, Ori used the technique of late time exp ansio n (ba¬ 
sically similar to the method presented in sec. VIII of the 
present paper), to ex plo re the late time behavior of scalar 
fields inside charged and rotating black holes. In 
this analysis, boundary conditions for the wave evolution 
were assumed on the event horizon (in the form of an 
inverse power-law in v), and the asymptotic late time 
{t ^ M) behavior of the wave was deduced inside the 
black hole, up to the inner horizon. This provided a tool 
for exploring the nature of the inner horizon singular¬ 
ity. With the results of the external analysis (generalized 
to charged and rotating black holes), a connection may 
be established between the form of the wave at the in¬ 
ner horizon to its form at null infinity, which, in turn, 
using the iterative scheme, can be derived as explicitly 
related to the form of arbitrary initial data outside the 
black hole. That would allow one, given initial data out¬ 
side the black hole (at large distance), to deduce the late 
time form of the wave at the inner horizon (including 
its accurate amplitude coefficient) without any assertion 
about the boundary conditions. 

One may also think of a more rigorous and coherent 
scheme, which includes the simultaneous analysis of both 
internal and external perturbations, in the framework of 
a generalized late time expansion. This generalization 
becomes natural when applying an expansion of the form 
(1^), as the coordinates v and r are both regular through 
the event horizon. Then a full treatment of both internal 
and external evolution is possible by fo llowing basically 
the same steps as described in sec. VIII, this time allow¬ 
ing the r coordinate to take its full range of values. 
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